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We study the effects of the equation of state on the nuclear suppression of heavy flavours in quark 
gluon plasma and estimate the initial entropy density of the system produced at the highest RHIC 
energy. For this purpose we have used the experimental data on the charged particle multiplicity 
and the nuclear suppression of single electron spectra originating from the semi-leptonic decays of 
open charm and beauty mesons. We have used inputs from lattice QCD to minimize the model 
dependence of the results. We obtain the value of the initial entropy density which varies from 
20 to 59 /fm'^ depending on the value of the velocity of sound that one uses for the analysis. Our 
investigation leads to a conservative value of the initial entropy density ~ 20/fm^ with corresponding 
initial temperature ~ 210 MeV well above the value of the transition temperature predicted by lattice 
QCD. 
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A thermalized system of quarks and gluons, called 
quark gluon plasma (QGP) is expected to be formed 
in the collisions of two nuclei at ultra-relativistic ener- 
gies [l| . Rigorous experimental and theoretical efforts are 
on to create and characterize this novel, deconfined phase 
of quarks and gluons. Lattice QCD (LQCD) calculations 
indicate that at a temperature ~ 175 MeV the entropy 
density (s) of the hadronic matter rises significantly due 
to the release of colour degrees of freedom which are con- 
fined within the hadrons at zero temperature. Therefore, 
it is of foremost importance to determine the value of the 
initial entropy density (si) / initial temperature (Ti) for 
the system formed in nuclear collisions at Relativistic 
Heavy Ion Collider (RHIC) and Large Hadron Collider 
(LHC) and assess whether the system is formed in colour 
deconfined phase or not. The focus of the present study 
is to estimate Si or Ti of the system formed at Au-f Au 
collisions at ^snn = 200 GeV. For this purpose we strict 
to take inputs from experimental data and LQCD calcu- 
lations to minimize the model dependence of the outcome 
of the present analysis. 

One of the possible way to estimate the value of the ini- 
tial entropy density is the extrapolation of the measured 
(final) observables backward in time through a suitable 
dynamical model. In absence of viscous loss the time 
reversal symmetry of the system is valid, therefore, the 
measured multiplicity at the freeze-out of the system can 
be used to estimate Si. The Si and the thermalization 
time (tj) are constrained by the measured (final) hadron 
multiplicity {dN/dy) by the following relation [2|: 



where Ax_ is the transverse area of the system can be de- 
termined from the collision geometry and k is a known 
constant (=3.7 for massless bosons like pions). The value 
of dN/ dy which is connected to Si through Eq. [T] is read- 
ily available for different collision centralities [i| . In Eq. [1] 
there are two unknown quantities, and Si both of which 



can not be determined from a single equation involving 
a single measured the dN/dy. Therefore, we choose an- 
other experimentally measured quantity the nuclear sup- 
pression of heavy quarks (HQ), Raa 3, which is sen- 
sitive to the initial condition and hence is very useful to 
estimate Si. 

The advantage of choosing the heavy flavour are (a) 
they are produced in early hard collisions and hence can 
witness the entire evolution of the QGP and (b) HQ 
are Boltzmann suppressed at the temperature range ex- 
pected to be achieved in heavy ion collisions at RHIC, 
therefore, the HQ do not determine the bulk features of 
the QGP. The magnitude of Raa depends on the amount 
of drag the heavy quark faces during its propagation 
through QGP. The heavy quarks being the witness of the 
early condition and the drag (and diffusion) coefficient is 
a temperature dependent quantity that makes the Raa a- 
good probe for the measurement of initial temperature. 

Earlier various attempts have been made to explain 
the experimental results on Raa - some of these are ad- 
dition of non-perturbative effects from the quasi-hadronic 
bound state Q , three-body scattering process @ , dissoci- 
ation of open heavy flavoured mesons by the thermal par- 
tons and inclusion of temperature dependent strong 
coupling Q. 

We briefly outline the procedure of evaluating the Raa 
for single electrons originating from the semi-leptonic de- 
cays of heavy mesons produced from the fragmentation 
of the HQ. The HQ while propagating through the QGP 
dissipates energy in the medium and hence its momen- 
tum gets attenuated. The magnitude of the momentum 
degradation gets reflected in the experimentally mea- 
sured quantity, Raa mentioned above. Theoretically the 
momentum evolution of the HQ in the expanding QGP 
background can be described by using Fokker Planck 
Equation (FPE) [9l-[l9|. The evolution of the probe i.e. 
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the HQ is described by the FPE: 
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7 and D are called the drag and diffusion coefficients, 
contain the interaction of the probe with the medium. 
To solve the FPE we need to supply the drag and diffu- 
sion coefficients and the initial HQ (charm and beauty) 
momentum distributions. 

There are two main processes through which the HQ 
dissipates energy in the QGP: (i) energy dissipation can 
take place due to the elastic collisions of the HQ with 
the quarks, antiquarks and gluons in the thermal bath, 
(ii) the radiative process due to which the HQ emits soft 
gluons (which subsequently get absorbed in the QGP) 
due to its interaction with the QGP. The details solu- 
tion of FPE with temperature and momentum depen- 
dent drag and diffusion coefficients including both the 
processes (i) and (ii) have been discussed in our earlier 
works [isj . The initial momentum distributions of charm 
and bottom quarks at RHIC energy {Js-^n = 200 GeV) 
have been taken from MNR code [20] • The method of 
solving the FPE numerically with temperature and mo- 
mentum dependent transport coefficients including other 
issues are discussed in [l^, therefore we do not repeat 
those details here. 

Classically the induced radiation takes place due to 
the jiggling motion of the propagating particle in the 
medium. Since the heavier particle jiggle less conse- 
quently induced energy loss is expected to be smaller 
[dead cone effect [2l| (see also [22, llSl)] for HQ com- 
pared to light particles. However, the experimental data 
from RHIC indicates similar amount of energy loss by 
heavy quarks and light partons in the measured kine- 
matic range. Various reasons like the anomalous mass 
dependence of the radiative process due to the finite size 
of the QGP [l^], development of dead cone due to high 
virtuality of the partons resulting from the dismantling 
of colour fields during the initial hard collisions 25| have 
been proposed as reasons for this observation. The au- 
thors in [26| concluded that the reduction in the energy 
loss of HQ due to radiative process is due to the dead 
cone effect but it is fair to mention that the issue is yet 
to be settled. 

The other effect which influences the radiative loss 
is the Landau-Pomeranchuk-Migdal(LPM) effect. This 
effect originate due to the interplay between two time 
scales of the system [23|: the formation time {rp) and 
the mean scattering time (t^c) of the gluon radiated from 



the HQ. LPM effect imposes certain constrain on the 
phase space of the emitted gluon [H, [l^ . Both the dead 
cone and the LPM effects have been taken in to account 
in evaluating the drag and diffusion coefficients in the 
present work (see for details). 

Now we discuss the QGP background with which the 
HQ interacts. The equation of state (EoS) plays a cru- 
cial role in describing the space time evolution of the ex- 
panding QGP from the initial state to the quark-hadron 
transition point. We use boost invariant hydrodynamic 
model [30] with the LQCD calculation EoS [31,] for the 
space time description of the matter. The velocity of 
sound (cs) as obtained in [3l| from LQCD calculations 
shows a significant variation with temperature (Fig. [1} . 
It starts with a very low value of at T ~ and then 
increases with T to reach the maximum value (c^ = 1 /3) 
corresponding to the value of a massless ideal gas. The 
EoS for almost baryon free QGP expected at RHIC en- 
ergy is take as: P = cj.e. The EoS sets the expansion time 
scale for the system as Toajp ^ [(l/e)de/dr]^^ ^ t/(1 + c^) 
indicating the fact that lower value of Cs makes the expan- 
sion time scale longer i.e. the rate of expansion slower. 
Therefore, for given values of Ti and Tc the life time of 
the QGP will be longer for smaller Cs . The value of Tc is 
fixed at 175 MeV. 

We solve the FPE numerically [3] for momentum de- 
pendent drag and diffusion coefficients to get the charm 
and beauty quarks momentum distribution. The solution 
then convoluted with the fragmentation functions [3^ to 
obtain the transverse momentum distribution of the D 
and B mesons which subsequently decay to create lep- 

tons [H,!!!]. 

In the same way the lepton spectra from the heavy 
flavours produced in p-p collisions can be calculated from 
the charm and beauty quark distributions which enter 
as initial conditions to the FPE. The solution of FPE 
contains the effects of drag (quenching) on the HQ where 
as the the initial distributions of HQ does not contain any 
such effects, therefore the ratio of these two quantities, 
the nuclear suppression factor, Raa act as a marker for 
the medium. This is observed experimentally through the 
depletion of Raa at high transverse momentum (pt)- 

As discussed above for larger Cs the expansion time 
scale is shorter i.e. the QGP life time is smaller, con- 
sequently the HQ spends less time in the QGP which 
ultimately leads to less suppression of the single electron 
spectra originating from the decays of HQ. Therefore, we 
take the following strategy to obtain the allowed range of 
initial temperature. We take the highest possible value 
of c^(= 1/3) for the space-time description of the flow- 
ing QGP background, in the present approach this will 
lead to the maximum value of Ti. In this case the HQ will 
spend the lesser amount of time in the QGP. Therefore, to 
achieve the experimentally measured Raa one will need 
larger drag or in other word larger initial temperature. 
The results for = 1/3 is displayed in Fig. [5] The value 
of Ti obtained from the analysis for this case is 300 MeV, 
the corresponding value of Si = 27r'^ g^f fT^ /AB ~ 59/ fm'^. 
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The value of Qef / 38 is extracted from the variation of 
s/T^ with T provided by the LQCD calculations [31)]. 

In Fig.[3]results for = 1/4 are depicted. For = 1/4 
the HQ spend longer time in QGP than for = 1/3. 
Therefore, with a lower initial temperature, Ti = 240 
MeV and Si ^ 29.34/fm^, the data can be reproduced. 
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FIG. 2: (colour online) Variation of Raa withpr for = 1/3 
and Ti = 300 MeV. 
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FIG. 1: Velocity of sound squared as a function tempera- 
ture fsl 



For Cg = 1/5 the HQ spend longer time in QGP (com- 
pared to the case when = 1/4). Therefore, the ob- 
served suppression dictates to reduce the initial tempera- 
ture. In this case the data is well reproduce with Ti = 210 
MeV and s; - 19.66/fm3. 

Further lowering of will make the value of Ti large 
(for given dN/ dy) enough to contradict other results like 
the observation of larg e hadronic elliptic flow which re- 
quires small Ti (see !355 for review). That will also results 
in lower Ti with which it will be difficult explain other 
experimental results. 

For all the theoretical results displayed in Figs. [2l [3] 
and |3]we have kept the quantity dN/dy constant conse- 
quently the value of changes to 0.6, 1.17 and 1.7 fm/c 
for Ti = 300, 240 and 210 MeV respectively The changes 
in Ti is forced by the change in the EoS. In Fig. [5] we 
show the variation of Ti with obtained by constraints 
imposed by the experimental data on Raa and dN/dy. 
The value of Ti varies from 210 to 300 MeV. In this con- 
text we compare the value of Ti obtained in the work with 
some of those reported earlier. In Refs. Q the value of 
Ti is obtained as ^ 375 MeV from the study of heavy 
quark suppression. From the simultaneous analysis of 
light and heavy quarks suppressions in Ref. [s^ a value 
of T, = 400 MeV is obtained. The authors in Ref. 
mentioned the values of the initial gluon rapidity distri- 
bution, dNg/dy, which may be converted to Ti — 290, 270 
and 310 respectively. 

It is interesting to note that the lowest value of Ti 
obtained from the present analysis is well above the 
quark-hadron phase transition temperature, indicating 
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FIG. 3: (colour online) Variation of Raa with pr for for 
1/4 and Ti = 240 MeV. 



the fact that the system formed in Au-|-Au collisions at 
y/syiN = 200 GeV might be formed in the partonic phase. 

In summary, we have studied the effects of the EoS 
on the suppression of single electrons originating from 
the decays of heavy flavours produced in Au-|-Au colli- 
sions at Y^SNjv = 200 GeV. We found that the initial 
temperature may vary from 210 to 300 MeV depending 
on the velocity of sound, which sets the scale for the ex- 
pansion that one uses. We have used experimental data 
(charged particle multiplicity and Raa of heavy flavours) 
and LQCD results (c^, gcff etc.) to keep the model de- 
pendence minimum. The effects of transverse expansion 
is neglected here. With the transverse expansion the HQ 
will (1) travel longer path (2) with diluted density. How- 
ever, the two competing effects (1) and (2) will have some 
sort of cancellation due to which our final conclusion may 
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FIG. 5: The variation of Ti with for fixed dN/dy. 



FIG. 4: (colour onUne) Variation oi Raa withpr for = 1/5 
and Ti = 210 MeV. 



not get altered. 
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